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In this Letter we study the entanglement of a quantum radiation field interacting with a dielectric 
medium. In particular, we describe the quantum mixed state of a field interacting with a dielec- 
tric through plasma and Drude models, and show that these generate very different entanglement 
behavior, as manifested in the entanglement entropy of the field. We also present a formula for a 
"Casimir" entanglement entropy, i.e. the distance dependence of the field entropy. Finally, we study 
a toy model of the interaction between two plates. In this model, the field entanglement entropy 
is divergent, however, as in the Casimir effect, its distance-dependent part is finite and the field 
matter entanglement is reduced when the objects are far. 



The theory describing electromagnetic fluctuations in a 
linear material medium is an essential tool in understand- 
ing matter-radiation interaction. Classically, thermally 
excited fluctuations as well as purely quantum, zero tem- 
perature, fluctuations in macroscopic media have been 
widely studied in terms of susceptibilities (see e.g. [UH]). 
While the use of such response functions have been highly 
successful in describing various phenomena, the actual 
quantum state of the field is usually not described. In- 
deed, even at zero temperature, a field interacting with 
realistic materials (as opposed to idealized Dirichlct or 
Neumann boundaries), is in general not in a pure state 
and should be described in terms of a mixed state density 
matrix. In modern terms, this mixed state is viewed as 
a consequence of matter-field entanglement, and carries 
a nonvanishing von-Neuman entropy, sometimes called 
"entanglement entropy" (EE). 

In this Letter we describe the quantum state associ- 
ated with a field in a medium given a dielectric function 
e. The effective action for such a field includes frequency 
dependence. Usually, such frequency dependent action 
is associated with response functions, however, here, we 
have to understand the action as specifying an instanta- 
neous form for the photonic density matrix. We do this 
by describing the effective density matrix of the field. 

As a measure of the field-matter entanglement, we use 
the entanglement entropy (EE) of the field. EE has been 
the subject of intense investigations in recent years, in 
a large variety of quantum systems, with applications to 
quantum information, condensed matter theory and high 
energy theory (For a review see e.g. [3]). In particular, 
the entropy of radiation coupled to matter has been con- 
sidered in numerous works. Usually the focus is on a sin- 
gle degree of freedom coupled to a bath of oscillators: For 
example, the entropy of a spin in the spin-boson model 
within the frame work of the Caldeira-Legget model [4 
was considered in [5] while the entanglement of a single 
radiation mode with an array of spins was studied in [5] 
for the Dicke model . On a different note, the EE of spa- 
tially separated intervals of vacuum (or ground state of a 
spin chains) has been considered in [8] . Here we con- 
sider a situation distinct from these works, and more akin 



to the scenario considered in studies of Casimir and Van 
der Waals interactions between macroscopic bodies: that 
of a field in contact with macroscopic dispersive bodies. 

Consider the electromagnetic field interacting with a 
dielectric medium at zero temperature. We assume that 
no external charges or currents are present and work in 
the gauge Aq = 0. The long wave-length effective action, 

S= ^ K j d 3 xdwA*(x)[cA(w,x) - V x V]A u (x), (1) 

encodes the field interaction with the material through 
the dielectric response function e(w,x). As we show be- 
low, the appearance of frequency dispersion in e(w,x) 
means that an action like ([I]) does not describe the state 
of the field as a ground state of a field Hamiltonian; in- 
deed, the action is nonlocal in time and so cannot be 
quantized to yield a proper quantum Hamiltonian. 

Clearly, one possible effect of coupling to the environ- 
ment may be thermalization of our field. It is therefore 
tempting to try and describe the instantaneous state of 
such a field as effectively thermal i.e. p ~ Z~ 1 e~ l3Hc tf for 
a reasonable effective field Hamiltonian H e f / . Such effec- 
tive "Entanglement Hamiltonians" have been the focus 
on recent studies due to their relation to the conformal 
edge spectra of fractional quantum hall systems [91 HH] ■ 

We show, however, that in the field-matter system, the 
field cannot be viewed as simply thermal. Rather, differ- 
ent ranges of momenta feel different effective tempera- 
tures. While the occupation probability of UV photons 
vanishes with high momenta, we still find that the field 
in ([I]) carries a logarithmic UV divergent quantum "zero 
point" entropy and show it's cutoff dependence in 

Finally, we consider, a la Casimir, the distance de- 
pendence of the entropy associated with field interaction 
with a pair of objects separated by vacuum, and show a 
toy model in which the "Casimir" entanglement entropy 
is UV finite and decays with distance. Physically, this 
means that the field becomes more entangled with the 
plates as they come closer. 

Here we study a simplified scalar field version of the 
electromagnetic field action 0. 

S = _L f d 3 .Tdc^C(x)[w 2 e(^,x) _ V 2 ]^(x). (2) 
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Note that the UV dependence study bellow is valid for the 
full electromagnetic problem, since in homogenous sys- 
tems separates naturally into two scalar fields describ- 
ing the two polarization modes. When the permittivity is 
independent of u, the action is local in time, and one can 
easily quantize the associated scalar action assuming the 
conjugate momentum ir<p can be expressed in terms of <p 
and doesn't depend on external fields. Such an action fol- 
lows from the Hamiltonian H = ^ j d 3 x[-^ + (Vc/>) 2 ]- 
It describes, at zero temperature, a pure state, and as 
such will have no entropy [32] ■ 

The situation is fundamentally different if e is uj de- 
pendent. The non- locality of the action ^ in time sig- 
nals that our system is coupled to external degrees of 
freedom which have been integrated out, yielding a non 
trivial temporal response kernel. In such a case, the sys- 
tem cannot be in a pure state implying that our radiative 
system is entangled with the matter fields. 

To proceed, let us briefly review the method of cal- 
culating entropies of Gaussian states (For details see, 
e.g. [HIE])- The calculation is facilitated by the fact 
that all information lies in the two point functions of 
the field. For a scalar field with n degrees of freedom 
cj> n and conjugate momenta n n , one defines the vector: 
(Oi,..C>2„) = (0i, 7ri, ..<j) n , 7r„). The state of the field is 
then determined by the covariance matrix: 



ljk =2Re((O j -(O j ))(O k -(O k })). 



(3) 



7 can be brought into a Williamson normal form 
W'yW 7 ' = diag(fj,i, ...fj, n , fii, ...fi n ) by means of a sym- 
plectic transformation W £ Sp(2n) preserving the canon- 
ical commutation relations [Oj,O k ] = i&j.k, where a is 
the 2n x 2n matrix 



o 1 

-1 o 



W f aW = a. 



The occupation probability of the normal modes are 
directly related to the symplectic eigenvalues //$ of 7. 
Finally, the density matrix can be written as: 



^-i e -Ei«i*i *i 



log 



f 



Pi - 1 



(4) 



where is a bosonic creation operators for the normal 
mode i. The entropy is given by: 



s = T,i hQii) ; h(/j.) 



2 



(5) 



An additional, technical simplification occurs, if we as- 
sume that no (</>7r) correlations are present. It is then 
known that the symplectic eigenvalues are square roots 
of eigenvalues of T = QT-L, where Q and % are field and 
field momentum two point functions, respectively. 

For the translationally invariant medium the sym- 
plectic eigenvalues can be labeled by momentum and 
given by p, k — 2n~ 1 (gkhk) 1 / 2 , where t/k^k is the 



product QV. in k space given by g k — (c/> 2 )k = 
/ dxe lx ' k (0(x,t')(?!)(O,t))| t /^ t+ and similarly h k = (ir 2 ) k . 
The field entropy per unit volume can be written as: 



S field = -V lr Trpf ie i d logp f ie i d 



dk/i(/x k ) (6) 



For concreteness, let us choose a typical dielectric func- 
tion e to use in ^ and We take: e(uj, x) = 
1 + 47rx(o;,x). With a typical susceptibiliy of the form 

(for a conductor we will add a Drude 
The action ^ allows us 



Xb 



function Xc = Xb + i ^L^ ] 
to compute, per k mode: 



gk 



duj- 



1 



w 2 e(k, iui) + k 2 
and, using time point splitting, 
1 



(7) 



,„ d ^Mk^) + P (fc 2 + (8) 

Using, the integrals ([7]),([8]) and Q we represent the den- 
sity matrix as p ~ exp(— X) k Wk< ^k ^k)- 

One may try to interpret this state as thermal, i.e.: 
p ~ exp(—/3H e f f). The effective field Hamiltonian is 
then given by: H e ff = J dxdx'u(x— x')$+$ X ' where u = 
[3~ x /dkitke jk '( x ~ x ). However, we find that generically 
the fourier transformed u gives us a non-local H e ff- 

Alternatively, we may interpret itk as Uk = fi k Ekj ree 
where E k j ree are the photon energies without the inter- 
action. Thus, different momenta k feel different effective 
temperatures. Estimating Uk from the integrals (fTb 5 (IsJ) , 
we find for the soft modes, Mk = /3 k E k j ree cx Vk and so 
the effective temperature: 



Vk 



CIS 



0. 



The energy and number of occupied soft modes per unit 



volume up to a given k m are proportional to k 



d+1/2 



and 



km 1//2 , respectively, are finite and small. The number 
variance of modes up to k m is (5N 2 ) ~ J d d kn k (l+n k ) ~ 
in d > 1. However, for d — 1 we find a curious infra- 
red divergence: (8N 2 ) = — log £/, where £/ is an infra-red 
cutoff, inversely proportional to the system size. 

Next, we find that the field entropy Q suffers from 
a UV divergence. To do so, we estimate the symplectic 
eigenvalues /ik to lowest order in lo p ,u c and use ((6J. We 
find that /z k ~ 1 + ■K~ l 2Wp r y p log k + .. for k 3> 1 , and 
integration over momenta yields (in 3d, h = c = 1 units) : 

(w 2 7p + w 2 7c ) log(A) 3 
S field = [ w 2 w log(A) 2 7c = 7 P = , (9) 

plasma model 

where A is a high momentum (UV) cutoff. Numerically, 
the approximations used in ^ actually recover the cor- 
rect A dependence even for large values of w p ,u; c , since 
the dielectric response decays at the large fc, oj limit. 
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It is interesting to observe the special place of the "pure 
plasma" limit response function. We can easily under- 
stand the result (|9| as follows. An idealized plasma sim- 
ply adds a finite mass to the field inside the region it 
occupies, completely expelling frequencies smaller than 
Lo p . Explicitly, substituting the plasma permittivity limit 

form x — ~ zfi m t ne action ^ , it produces a mass term 
for <j). Thus, the resulting action is consistent with a 
Hermitian field Hamiltonian, and as such, at zero tem- 
perature, to a pure state. Interestingly, the use of a pure 
plasma in the computations of Casimir energy and En- 
tropy has been at the heart of a recent debate [T3Tll6| . 
note that the distinction between the Casimir entropy in 
the two models is manifested in the full quantum entropy 
computed herein [33] - 

Having shown that the entropy Sfi e id is UV divergent, 
it is natural to ask, in analogy with the Casimir effect, 
what is the distance dependence of the entropy of inter- 
action with two distinct bodies A and B1 Is it UV finite? 
To answer such questions, we define: 



S R (A, B) = S(A U B) — S(A) - S(B), 



(10) 



Sr should not be confused with the Casimir entropy, de- 
fined as Sc {A, B) = — lim T _j. 9tF c , where F c is the 
Casimir free energy, obtained by subtracting all distance 
independent terms from the free energy of the EM field 
in the presence of the bodies A, B. Also, Sr is distinct 
from the "relative entropy" of probability theory, as we 
are comparing different systems, and not merely different 
statistical information about the same system. 

The relevance of Casimir entropy Sc to understand- 
ing thermal corrections of the Lifshitz formula has been 
pointed out in many papers (see e.g. [13l [T71 [18] ) , where 
it was noticed that as T — > 0, Sc may not go to zero 
when using the Drude model (as might be expected by 
the Nernts theorem). It is interesting to note that while 
the Casimir EE Sr is distinct from Sc, a similar behav- 
ior is observed in (j9J). In addition, at high temperatures 
we expect Sr = Sc, as most of the field entropy will be 
thermal (Technically, the relevant Green's function gets 
its major contribution from the oj — Matsubara pole), 
in addition, the Casimir force is entirely entropic |19j . 

In recent years, numerous results for Casimir interac- 
tion between objects have been obtained using a TGTG 
representation of the energy [20]. In these methods 
one separates between T operators associated with lo- 
cal properties of each body, and a free green's functions 
interpolating between them. See, e.g. ([HH23]). A nat- 
ural question is: Can we describe Sr in similar terms? 
Using the methods of [23], assuming that <j)ir correlators 
vanish, we find an analogous expression for Sr (A, B): 



sr(a,b) = -±j; 



X/2+i. 



dx- 



log 



I] 



Trlog 



( 



/2-ioo "**" 

1 - T^KI ( K aKb + K AlJB -K A - K B ) 



£(0J) 



£ = 1 



E(OJ) 



FIG. 1: "Casimir" entanglement entropy Siz of plates. 



Here K A = x _ 1/i+is {T A - 1/4) is defined via T(x, x') = 
i J((f>(x)(j)(x"))('K(x")Tr(x'))dx" , and similar expressions 
hold for Kr, Kaub- Expression (111 is similar to the 



TGTG formulas in it's form: an integral over the TrLog 
of a combination of Green's functions. It differs from such 
formulas in several aspects: The integration variable x is 
not a frequency variable, but rather an auxiliary spectral 
variable, the presence of the term Kaub doesn't allow 
for full separation into local object properties and free 
propagators and the non-analyticity at x = 1/4. All of 
these make the formula harder to use than the TGTG 
formulas. Nevertheless, it is a useful starting point for 
multiple scattering expansions. 



Here, we defer a detailed study of the formula ( 11 ) and 



proceed instead to study the dependence of entanglement 
on the distance between bodies in a simple toy model: 

Consider a field interacting with a medium, defined 
through the following (Wick rotated) action: 



where tp is a matter field which is confined to the body(s) 
B. This action may be viewed as a purification of the 

v = T^|fi)<n| 

tjj space. The 

2 

= — — 2 of the 
response of to a transparent, but dispersive medium. 

Finally, to render the problem exactly solvable, we will 



state p0 of 4>, i.e. we write our state as 
for some pure state |fi) in the larger (j. 



action ( 12 ) corresponds to the form x(w) 



assume that the 4>ip coupling in ( 12 1 comes from the form 
uj p J tpcj) in the classical Lagrangian. Note, that while on 
the classical level this Lagrangian coupling is equivalent 
to — u! p J cfnp, the actual quantum mechanical transfor- 
mation required to affect this for gauge fields is of the 
Power-Zineau-Wooley (PZW) type used to interchange 
the minimal coupling A ■ J to E ■ P in quantum elec- 
trodynamics |24j . However, PZW mixes radiation and 
dipole degrees of freedom and does not preserve entan- 
glement properties. 

To compute the entropy of the field <j> more efficiently, 
we use that S^ = S^ (noting that ( 12 1 describes a pure 
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state). <S^ is obtained from the effective action for tp: 

S cS m = f B d 3 *dV[c^ 2 (x|G fl |x') + 

5(x-x')4^V]^*(x)^*(x'), (13) 

where Gq is the restriction of the free Green's function 



to the bodies B. 



~~ -w 2 +V 2 _ 

The two point function of tp is: 



(V(x)V(x')) = / <x| 



W)P (14) 



with a similar expression holding for the conjugate mo- 
mentum of tp (obtained from the Lagrangian in Spure): 



i 

' 167T 4 



5 _2 rUp w lj- 



71"V(x) 7r V(x')/ 

8^<5(x-x')|x')^ (15) 



Next, we study the entropy generated by parallel 
planes separated by a distance L, as illustrated in Fig(T] 
To control UV behavior we place the system on a lat- 
tice, and use the discrete analogue of the free Green's 
function (in ID) as: Go(n,m,k±) 



-q , 



2 i .2 

cosh q = 1 + - 2 -*- (k± is parallel to the plates). 

We first consider the ID case. In ID, the planes reduce 
to a pair of sites, and the kernel ( ^ — h uj p 2 uj 2 Gq)~ 1 
in eqs. ( 14|15 1 is a 2X2 matrix, amenable to an analytic 
treatment. Estimating the two point functions (14 1 and 



2 sinh q 



where: 



(151 using Watson's lemma, we find, at large distances 



L>1 (assuming loq > uj p and up to order i 



Qnm = (ip(ri)ip(m)) ~ A_ 



A 

32 7T 



32 iv 



G 

4, . 2 




-4 r 



_ 4^ _ 327T 4 i^(^ arccos[^]-2 v /^g-4) 



Wo 2 K 2 -4)3/2 



(16) 



Similarly, the momentum correlation function H nm is 



rLr. 



-l 

32tt 5 



192tt 4 ^ 
" 4 £ 4 



192ir 4 



c 4 L* 



(17) 



4tt 



32Tr 4 ^ (2y / ^ 2 -4+(wQ-8) arccos^]) 



The dominant terms in the product QTi. are proportional 
to = -^qAqAh, and are independent of distance. 

As expected on physical grounds, the entropy strongly 
depends on the pinning frequency ujq. We note that 
the entropy is maximized when loq is small as possi- 
ble, i.e. loq ~ Lu p . In this limit the tp field is only 
weakly restrained to tp = values, and thus generates 
large entropy for the <p field. In the weak coupling limit 
w p « 1 we find the eigenvalues of QH to behave as 
9 (ujq — LJ P ) ~ 4 + 27TCj p log In the opposite limit, 
when wo 3> w p the ?/> field is essentially "pinned" down, 



00 we 



and we have 9 ~ 
recover the minimal possible value of the product QT-L 
allowed by the uncertainty principle. 

Finally, the symplectic eigenvalues of the covariance 
matrix are found, to leading orders in j^,oj p , to be 

\± = Ve± AHUpi^oj^VeL 2 )- 1 (is) 

and the entropy using (|5| for large L to behave as: 

= S{L ->■ 00) + 32vr 7 a;2 a ; -3 L -4_ ( 19 ) 

At dimensions D > 1, the momentum is a good 
quantum number. We proceed to compute the en- 
tropy per k± by estimating (tp(x,k±)tp(x' J k±)} and 
( 7T ip(x,k ± ) 7r if)(x r ,fcj_)) ■ We find that the L dependent, off- 
diagonal elements in ( 16|17 ) are exponentially decaying 
as e _L ' fci '. Integrating over k^_ yields 



50 = S{L -> 00) + C4(cjo,lj p )L 



-3-d 



(20) 



where Cd(wo,^p) is a constant. 

Discussion: We studied the quantum state of a field 
in the presence of a dielectric material. We found that 
such a field is described by a density matrix whose Von- 
Neumann entropy diverges as described by eq. |9|. The 
state cannot be considered as thermal, but rather the 
photons have a k dependent effective temperature. The 
occupation number of modes behaves as ~ 1/Vk per 
unit volume for "soft" photons. A similar situation may 
arise when considering the phonons in solid. This effect 
is reminiscent of the infra red problem in quantum elec- 
trodynamics where, infinite numbers of soft photons are 
generated in transition amplitudes (see, e.g. [25 ). The 
situation here is different in that the system is in equilib- 
rium, however we are integrating over dipole transitions 
in the material. 

To find long distance features, unaffected by UV di- 
vergence, we considered Sn, the distance dependent part 
of the entropy of two objects. In eq. (11) we present 



a novel general formula for the "Casimir Entanglement 
Entropy" . Finally, we considered a toy model for the en- 
tropy generated in the field due to interaction between 
two thin planes. In this model S-n shows a L~( 3+d ) decay: 
Thus, the field matter EE is reduced when the objects 
are far. Further investigation is needed to determine the 
power law in more realistic models. 

The measurement of the full quantum entropy is, in 
general, quite hard, but possible to achieve in certain 
situations [25l429| . In the case of Gaussian fields with 
translational invariance, it can be computed from mode 
occupations. We mention that many of the properties 
discussed are also valid for other fields, such as phonons 
in a solid. For example, occupation numbers of phonons 
in ultracold atoms or trapped ion systems may be mea- 
surable in experiments [301 131] , and may reveal some of 
the features described here. 
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